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ABSTRACT 


In this report we analyze the transpiration cooling by He bubble injection of the long LOx 

tank feedline heated by the environment heat. We consider possible hazards that can arise in 

the proposed design of the SLS core stage where the feedline length is much longer than that 
used in the Space Shuttle. 
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Nomenclature 


p 

= density 

P 

= dynamical viscosity 

p 

= pressure 

Pr 

= Prandtl number, Pr=p.Cp/k 

Pc 

= critical pressure for gas 

Ra 

= Raleigh number 

T 

= temperature 

Re x 

= Reynolds number 

T c 

= critical temperature for gas 

Gr 

= Grashof number 

T s 

= surface or interface temperature 

Nu 

= Nusselt number 

u 

= velocity 



c 

= sound velocity 



M 

= Mach number, M = u/c, M 0 = u/c 0 

Subscripts: 

C v 

= specific heat for the constant volume 

g 

= gas 

Cp 

= specific heat for the constant 

L 

= liquid 

y 

pressure 

env 

= environment 

g 

= ration of specific heats; y= Cp/c v 

t 

= tank 

R 

= gravity 

w 

= tank wall 

R« 

= radius 

ull 

= ullage 

s« 

= radius of tank 

eva 

= evaporation 

Ssu 

= surface of tank cross-section 

P 

condensation 

r 

= total surface of tank 

con 


Ld 

= thermodiffusion length 

d 


Ql 

= heat of evaporation 



h T 

= heat transfer coefficient 



j 

= mass flow density [j = up] 

Abbreviations: 

J 

= mass flow [ J = upS] 



q 

= heat flow density 

LH2 

= Liquid Hydrogen 

Q 

= heat flow 

GH2 

= Gaseous Hydrogen 

Qc 

= convection heat flow 

MNM 

= Multinode Model 

Qr 

= radiation heat flow 

SLS 

= Space Launch System 

Pr 

= thermal expansion coefficient 

3D 

= Three dimensional 

5 

= boundary layer thickness 



k 

= thermal conductivity 



V 

= kinematic viscosity 




Introduction 

Here we study main physical processes taking place in LOx tank during the pre- 
pressurization stage and the rocket flight. The goal of this research is to analyze the risks 
associated with the functioning of the external LOx tank of the Space Shuttle (Fig. 1 ) during the 
pre-pressurization and in the course of the flight and their scaling to the proposed SLS design 
with 5 engines (the situation with 4 engines is less critical). The risk assessment for the 
replenish stage is presented in a separate report. 

We use both the protocol [1] and the telemetric data [2] for Shuttle flights to fit an uncertain 
parameter of the multinodes model (MNM) developed in our group using results of 
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H : =15.05m 


Fig.1. Design of the Shuttle tanks: the shape and parameters of the LOx tank 
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Fig. 2. The shapes of the LOx tanks analyzed in this work 

works [3-5]. After an appropriate rescaling the model is further applied to study the temperature 
stratification, the ullage pressure dynamics and the associated risks in the LOx tank in the 
framework of the proposed SLS designs. This constitutes the risk scaling problem. Two different 
designs of the LOx tank will be considered. The first is tank of the shape of LOx Shuttle tank 
(Fig.1), but with elongated central part to contain 5/3 of the Shuttle tank LOx volume (Figure 2, 
left). The second is tank of the shape of LH2 tank of the Shuttle, rescaled to contain 5/3 of the 
Shuttle tank LOx volume (Figure 2, right). 


1. Main physical processes and pathways of heat transfer in the LOx 

tank 

The general heat balance in the LOx tank can be written as 

Qinjecl Qenv Q ullage +&*,+&• 0 ) 

Here Q mjea +£? t ,„ v is the total heat added to the tank as a result of hot gas injection from the 
diffuser and the interaction with environment. This heat is spent on heating the gas in the ullage, 
Q u „ age , heating the tank wall, O wall , and on condensation-evaporation processes including 

heating of the surface layer of the liquid , Q, . The heat flows associated with these pathways 

are sketched in Fig. 1, which also defines the variables of the physical model. 

We now proceed to analyzing these pathways and estimating their relative importance for 
the problem of temperature stratification and pressure dynamics in the ullage. 
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Fig.1-1. Main physical processes, variables and pathways of heat transfer. 


1 . 1 . 


Condensation-evaporation processes 


The timescale of the ullage pressure changes is large (~1 sec) compared to the timescale of 
thermalization of the liquid-vapor interface, which can be estimated as follows [6]: 

= (27r 2 K L C L p L R vapor Tl ) / (pl{T, )q; ) < 0. lmsec (2) 


* thermo 


Her ep, g ,C I j; ,/c Lg are density, heat capacity, and thermo-conductivity of liquid (vapor); 
Rvapor gas constant; q L - heat of vaporization. As a consequence, the quasi-equilibrium near the 
interface can be assumed for all times with the surface temperature T s determined by the 
vapor pressure p vap0l/ r = Pgox [6,7]: 


T=T 


Pc, Ox'' 


\n 


p.cn=Pc 


f j\ 


(3) 


Here T c and p c are the critical temperature and pressure of the gas, p s is the pressure of the 
saturated vapour and the constant X=7.2 for oxygen [7], Note, that p s =1atm at T=T L =90K, during 
of the prepress and flight p GQx >1atm and T_. >T L (T s (p GOx =2.25atm=36psi)=100K > T l =90K). 

The mass flow at the interface is limited by the heat balance_between the hot gas and the cold 
liquid surface. In what follows we analyze the condensation-evaporation processes during the 
pre-pressurization. Condensation-evaporation during flight is discussed in Section 3. 

The condensation of GOx prevails over the evaporation during the first stage of pre- 
pressurization (Fig. 1-2, -152s < t < -138sec). This result can be understood as follows. In the 
beginning of He injection (t= -152sec, see Section 2) the Ox vapor is compressed toward the 
liquid-gas interface due to the relatively high density of GOx. As a consequence the partial 
pressure of GOx increases (Fig. 1-3). 
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Fig. 1-2. GOx condensation mass flow and total condensed mass. Results of MNM 

simulation. 

Accordingly, the condensed mass accumulates at this stage and reaches the maximal value 
of ~1kg which supports the estimations presented in the text. At the later stage (-138 sec <t<0 
sec) the condensation rate is seen to oscillate due to the pressure modulations and decreases 
on average as predicted. The average condensed mass decreases at this stage. 

The 3D simulations confirm this conclusion (Fig. 1-4). The growth of the GOx partial pressure 
stimulates condensation of the compressed GOx, accompanied by the rise of surface 
temperature 



Fig. 1-3: Partial pressure dynamics of GOx and GHe during the pre-pressurization in 
the ullage layer near the liquid-gas interface. Results of MNM simulation. 



He Flow 
T=316K 


Inflation 30 Layers 
First Layer - 0.1 mm 




Fig. 1-4. Stratification mass fraction of GOx and He during the pre-pressurization 
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in accordance with Eqs. (2). The liquid surface temperature T s is determined by the partial 
pressure p GOx of Ox vapor near the interface: T s >T L at p GOx »'1.75atm > p s (T)=1 atm, where 

T l is the liquid temperature far from the interface. The temperature profile near the interface 
and the heat flows are schematically shown on Fig. 1-5. 



Fig. 1-5: Temperature profile near the liquid-gas interface 

Since T s >T l and T >T S the convection of both the liquid and the vapor is absent. 

Moreover, our simulations have also shown that the convection induced by the hot gas injection 
from the diffuser and natural convection near the wall is indeed very small near the gas-liquid 
interface. Thus, the heat balance in condensation-evaporation processes is established mainly 
by means of thermo-diffusion for the given surface temperature T s ( p GOx ) . The solution of the 
one-dimensional heat equation 


dT _ d 2 T 
l ' sPl ' % dt ~ Kl * dx 2 ’ 

for the boundary condition T(0,t ) = T S = const, T(x, oo) = T L is given by 


(4) 


T{x,t)=T L +(T s -T L )erf 


* k lJ 


(5) 


Here p, ,Q g ,K, g are density, heat capacity, and thermo-conductivity of liquid (vapor); x=0 

corresponds to the liquid-vapor interface, and erf is the error function. Using Eq. (5) we can find 
the heat flux from the warm interface into the liquid: 


Q l (x,1) = -k l 


dT(x,t ) 


dx 


= ~K, 


T,~T, _ \C,p,K L (T,-T L Y 


7Ct 


L DL = t \^ £ . (6) 


.r=0 ^DL ' /1 ‘ \C lPl 

Here L dl is characteristic thermo-diffusion length in the liquid. This makes it possible to 
estimate the condensation mass flow j cond from the following balance equation: 


fiw - {q L * r,„ (T t -T,))= Q,_ -Q,= S, 

J cond(t') 


( 


r -T T -T 

" 1 s 1 L _ g l 


^ £flt(0 L Dg(0 

>]K]Cj pj (T s -T l )- K GO x C GQx p GQx ( T g —T s ^ 


(7) 


{c[ L + C GOx ( T g T s ^ 
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The first term in the left-hand side of Eq. (7) describes the condensation flow rate and the 
second one describes the evaporation flow rate. The total mass of the condensed vapor is 
given by 

I 

Mcondi 0 = jX {t)j cond (t)dt , M cond (t) oc 4t at S, (/) » const (8) 

o 


According to Eq. (7) condensation will be dominant whenever 

(T, -TXPoo,)) < (rAPoa.)-TjJ -^£i , r,(p a »)= t, 

V ^aox^'ooxPoox 


n Y' 7 

Pgox 
K Pc J 


(9) 


where for the estimation of the thermo-diffusion length in the gas we can use parameters of 
the pure Ox, which is prevalent in the vicinity of the interface as argued above. We note that 
Eqs. (7)-(9) are valid during the flight as well (see Section 3). 

At the early stage of the pre-pressurization, condition (9) is satisfied. As argued above, the 
partial pressure of the compressed GOx vapor is controlled by the pressure of the injected He, 
PGOx~1-75 a tmi during the early stage of pre-pressurization (Fig. 1-3). In this case the surface 
temperature TJp=1.75)=97.5 > T l =90K and the condition (9) holds provided T K <450 K, 

which is certainly satisfied (temperature of injected He is 320K). Estimation [Eq. (8)] of the total 
condensed mass during the early stage of pre-pressurization ( ~14sec ) gives AT conf/ «l kg. 


This estimate is supported by the results of the simulation (Fig. 1-2). Since the initial mass of 
the vapor is = Ar*K.iia ge ~17kg and the mass of injected He is about 6kg, the condensed 

mass constitutes about 4% of the total mass of the gas in the ullage. (We note that the mass 
fraction of the gas involved in the condensation-evaporation is considerably smaller during the 
flight, see Sec. 3). 

At the later stage of pre-pressurization (-138sec < t < 0) the vapor becomes mixed with hot 
He and the GOx partial pressure will drop (Fig. 1-3). The gas temperature near the interface 
Tg~100K (Fig. 1-6). Therefore, at this stage the condensed mass decrease. The foregoing 
theoretical predictions are verified by the results of MNM simulation plotted in Figs. 1-3 and 1-6. 



Fig. 1-6. Ullage and wall temperature during the pre-pressurization: upper red and blue 
curves for the layers near top and bottom of the tank, respectively. Results of the MNM 
simulation. 

We note that at the later stage of pre-pressurization weak modulations of the partial 
pressure of GOx (Figure 1-3) are observed in Fig. (1-3). These modulations are due to 
pulsations in the hot He injection from the diffuser (see Section 2). The pulsations have been 
neglected in our estimations but the average decrease of the partial LOx pressure in the vicinity 
of the interface supports the theoretical argument. The average decrease of the partial pressure 
is accompanied by the decrease of the surface temperature (Fig. 1-6) which displays 
modulations as well. 

The condensation-evaporation during the flight is discussed in Section 3. 
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1 . 2 . 


Collapsing of bubbles during the initial stage of the pre- 
pressurization 


Ullage pressure increases from latm to 2.45atm during the pre-pressurization (-152sec <t< - 
138sec) due to injection of GHe. Therefore, the LOx becomes strongly overcooled and the 
bubbles which can be present in the liquid start to collapse. The characteristic time of collapsing 
can be estimated from the following heat balance equation 


S bubble 


\e .bubble ) 7, c , 

^ bubble A .( )x C ! /, bubble ’ 


'bubble 


bubble (0 “ bubble , oO , )> K 


1 col laps 


collaps 


_ r bubble,oP.s,Qx^i 


( 10 ) 


Here T sbub is bubble surface temperature that is determined by the ullage pressure 


T = T 

1 sjbubble * c 


P ullage 

. Pc 


v 1/7.3 


’ ^ collapse = 0 - 2seC 


r r 

bubble, o 


l mm 


(11) 


t collapse ~ 1-25 (5) sec at r buhhle o — 2.5(5 )mm. 

Here we took into account that p u ii a ge = 2.45 atm during the pre-pressurization (-138sec<t<0). 
Then, T s . bU bbie =102K and the density of saturated vapour in the bubbles is p s0x = 3.48kg/ m 3 . 

Validity of estimations given by Eq. (10) was confirmed by numerical simulations based on a 
more elaborate model of bubbles’ collapsing in overcooied liquids [8], Duration of the pre-press 
is longer than lOOsec. Thus, all bubbles are expected to collapse by the start time t=0. 


1.3. Heating of liquid in the boundary layer near the interface 


The liquid surface layer is heated due to condensation. It can see from Fig. 1-3 the partial 
GOx pressure near the interface p G ox < 1-7atm. Therefore, the surface temperature of the liquid 
according to Eq. (3) is T s £ 97K (T S -T L < 0.7K). The width of this warm liquid layer near the gas- 
liquid interface is determined by the thermodiffusion length of the liquid which increases with 
time and approaches a value of L dl =6mm toward the end of the flight. 

The upper bound on the mass of heated liquid can be estimated as 
Mm.* ~ S(t Jllghl ) p L L di (tjbgb, ) ~ 400Ag- . This value is negligible compared to the total LOx mass in 

the filled tank, 627,500kg. 

The heat flux into the liquid associated with the condensation is given by 


& - d< j; - k,. 

dt L DL {t ) V 7rt 

The total heat associated with the condensation of GOx can be estimated as 


( 12 ) 


Ql( t) * 2(7; - T L )S L Jt K L C L p L ! x, Q L (t mght ) « 10 7 J 


(13) 
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Fig. 1-6: Heated layers near the interface and the tank wall 


The total heat of GOx injected during the flight (t night =500sec) is 

Qinjeci = C pOx^ind in/ flight = 3 . 4 - 10 J , 

D — CM T — C — ^ P ullage ( fhgh ) ^tank _ a c inS r .. .. 

ullage ~ ^ pOx JV1 GOX 1 ullage ~ ^ pOx D ~ . - Z..HJ MU J V '^1 

k gox y ~ 1 

J mj = 1 .3 kg / sec, T mj = 800/? = 444/C 

We conclude that the relative heat content of the heated liquid layer is a small: 

Ql (J flight ) / Qirfccl = 3 % ■ 


1 .4. Heating of liquid in the boundary layer near the wall 


Heating of LOx near the tank wall occurs due to the environmental heat 
flux, Q anir — ^.„v„^u.nk which causes free convection of the liquid near the wall. The convection 
flow is turbulent and can be described by the following relations, discussed in a separate report: 


Lw 


<t., - =h(T„-T L ) = K, Ztt-M . h„ = 0. 1 5 K, 

n _ &Pi£jPl ( T l „ ~ T l )x _ 

" a x ~ rr - h.w — 

M l K l 


1/3 


(15) 


0.15 (Ra x ) 


1/3 


« 0.5 mm 


where /, is effective thickness of the heated layer near the wall. Here we use correlations for 


the turbulent free convection flow as in this case Rayleigh number Ra x >3x1 0 12 for x>1m (x is 
the elevation of liquid layer counted from the tank bottom). We note that the values of l Lv and 
heat transfer coefficient h w do not depend on x. According to Eq. (15) the temperature 


increase of LOx layer near the wall is less than 0.7K : (T Liy -7}) <0.7 K for typical 


environmental heat fluxes, q emir =(l30-^200)fT/w 2 . The total mass of the slightly heated liquid 
in the boundary layer can be estimates as follows: 
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( 16 ) 


M lox,t* 9 o. 7 k ^ S^p L l Lw ,S tmk - tank surface area, 

M lox,t* 9 oik s 2 1 $kg (360kg) for 3 (5) engines. 

We note that the total LOx mass in the tank is 627,500kg (1,046,000kg). Thus, the mass 
fraction of the heated liquid (T lw &90.1K) near the wall is negligible quantity and does not play 

a role in the ullage temperature stratification and pressure dynamics. Moreover, the thin warm 
liquid layers with T < 90. 7K near the interface and the tank wall, i.e. a slight stratification of LOx, 
does not lead to breach of condition NPSP = p lJI +p head -p fiic -p svap (T s )> 30psi for the net 
positive suction pressure due to high hydrostatic head pressure P he ad (see below Sec. 5.1). 


1 .5. Heating of the tank wall by hot gas: main mechanism of the gas 

cooling 


The heat flux from the ullage at T g to the tank wall at T w can be written as 


V,. = 


dq 

*7 v 


0.15 


= — = h gw (T g -T w ),h gw =—K g (T g ){Ra x ) *20 


dt 


W_ 

m 1 K 


Ra _ 8fi x rC ll P g (P^XT ll -TJx 3 




(17) 


K„ 


l lw = * \mm for p = 2Aatm, T, = 425 K, T w = 250 K. 

K 

Here we use the correlation formulas for the turbulent regime corresponding to large 
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Rayleigh number Ra>7x 10 for x~1m. The total heat flow is 

4SL _ , 


-f- = M r « - w,(0 = 


(18) 


The temperature of the ullage and the tank wall depends on the coordinate x along the wall. 
Temperature dependences of the specific heat C w for Al alloys used for the tank is shown in 
Fig. 1-7. 



Fig. 1-7: Specific heat for Al Alloys 2195 (left) and 2219 (right) 

These dependences were used for our MNM numerical simulations. For the estimation we 
assumedr =const, C W (T W ) = (C w ) = const and h = const . Then Eq. (18) implies 


ll 


( 19 ) 


dT w (T-TJ 


dt 


r gw(T w ) 

_ (C w )p w d„ 


gw 


= T g -(T g -T w0 ) Q xp(-~~), 

L gw 

■,r^»150sec at { C w ) » 600J/kgK. 


The changes in the temperature of ullage gas and in the temperature of the tank wall are 
shown in Fig. 1-8. The solid lines are plots of MNM-simulated time-traces of temperatures near 
the top (red) and bottom of the tank (green). The dashed lines are MNM-simulated time traces 
of the tank wall temperature near the top (red) and the gas-liquid interface of the tank (green). 
The blue solid lines give the estimated tank wall temperature near the tank top and the interface 
according to Eq.(13) for T gw =150sec. These analytical curves correspond fairly well to the results 
of the MNM simulation (Fig. 1-8). 



500 


g 400 
=> 

2 300 

CD 

Cl 

i 20 ° 


100 


40 60 
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Fig. 1-8 Time traces of the ullage and wall temperature: solid curves are MNM 
simulation results and dashed curves are analytical results 


It it follows from Eqs. (12), (13) and (14) the heat transferred from the hot gas to the tank wall 
is 

Qwall = Qinject + Qenv ~ Qulhge ~ Ql ~ 6.8 • 1 0 J , 

a„„ = CV.VVfc*. = 3.4 - 10 V, ( 20 ) 

&*, - =2.45-1 0 V, a.. = * 3.7- 10V 

Thus we conclude that heating of the wall is the main mechanism of gas cooling in the 
ullage. 

2. Pre-pressurization regime and fitting the Shuttle data 


The values of main parameters at the end of replenish (EOR) and at the engine start 
command (EST) are shown in Table 1 (from [1]). The variation of ullage pressure during the 
period, t = -300sec to t = 0, is shown in Fig. 2-1 [1], The nominal parameters of liquid and gases 
supply for STS -115 are presented in Table 2 [1], The replenish stage ends at t = -300sec, i.e., 

300sec before the start. The pre-pressurization starts at t^ press = -152.8sec with injection of 
the hot He. The nominal He flow rate is r Hi ,=0.3lb/sec = 0.136kg/sec and He temperature is T= 
572.4R=318K (see Table 2). 
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After the EOR one observes a leak of LOx which increases abruptly and reaches the rate of 
18lb/s = 9.1 kg/s (Fig. 2-2). Therefore, the ullage volume at can be estimated as follows: 

y( 0 ) = yEOR +8 1kg/ s A/ / p L = 5m’, A/ = (300- 152. 8> = 157.25 

where the liquid level =1.09 m is measured from the top of the tank downward. Hot He 

is supplied during the time A =(152.8-138.8)sec = 14sec until the pressure reaches 

20.54 psig=2.42atm (Fig. 2-1). At ^^=-138.8560 (liquid level xffl = 1.406m) the ullage 
volume 



Time, sec 


Fig. 2-1: Time-trace of the ullage pressure (STS 115). 

increases toy u ( j> =V u f c -9— 138.8sec/p ; =8w 3 - We also used the fact that the ullage volume 

sec 

increases with the pressure approximately as v / (p) = v i ( j >) — -^—(p-Po) due to the tank 

P\ ~ Po 

wall elasticity and the lowering of LOx level. 

The temperature, partial pressure and densities of GOx and He stratify at 
mispress - -138.8 sec as shown in Figs. 1-3 and 1-4 (MNM simulation) and Figs. 2-3 and 2-4 (3D 
simulation). 


Bleed flow of LOx from ET 



Fig. 2-2: Time-trace of leakage from LOx tank during stages EOR through ESC (STS 

115). 
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Table 1: EOR through ESC (from t=-300sec to t=0) including pre-press (STS-115) 


EOR 



ESC Through TO 


Ambieur Pressure: 


14.62 psia t 

L02 Losses Prior to ESC: 

5,904 lb*, 

Nose Cone AP 


0.01 psi 

GHe (Anti-Ice) Mass Added. 

1 lb m 

Ullage Pressure: (above ref) 


0.64 psid 

GHe (Prepress) Mass Added: 

131b* 

Ullage Pressure: 


0.65 psig 

Ullage Pressure: 

22.0 psig 

Ullage Pressure: 


15.27 psia 

Ullage Pressure: 

36.62 psia 

Liquid Level: 


412.58 Xt 

Ullage Volume « ESC (X=1.514m)321 ft 3 9.09m 3 

Ullage Volume: (X=0.929m) 3.879m 3 

137 ft 3 

Ullage Mass @ ESC 

49 lb* 

Ullage Mass: 


36 lb* 

ET Liquid Volume @ ESC 

19.354 ft 3 

Total Tank Volume 

555.51m 3 

19.617 ft 3 

ET Liquid Mass @ ESC 

1.377.432 lb m 

Bubble Volume: 

1.133m 3 

40 ft 3 

Press'd L02 Density 

71.17 lbjft 3 

Liquid Volume: 

550.5m 3 

19.440 ft 3 

Predicted L02 Usage to T-0: 

9,530 lb* 

Vented L02 Density: 


71.16 lb*'ft 3 

L02 Transferred to SSMEs: 

1721b* 

L02 Propellant Load: 

627.5 T 

1.383.3361b* 

Loaded L02 at T-0: 

1.367.730 ib ra 




Predicted Ullage Volume - % of 

1.66 

XT375 — ET top 



Liquid Voiume 


XT 3 9 6 .7- 0 . 5 m down 



Required Ullage Volume - % 

1.59 

XT4 12 - 0.95m down (initial liquid level 

= ullage) 

of Liquid Volume 



The ullage pressure drops at / > t ( p ^ press >-138sec due to cooling of the hot injected He (see 

Sec. 1.4) and expansion of the ullage volume. Pulses of He injection are applied when the 
ullage pressure drops below the threshold value of p th =20.54 psig=2.42atm (nominal value, 
Table 2). The results for the pressure time-trace from 3D simulation are presented in Fig. 2-5. 
Rectangular He pulses of nominal duration 0.5sec and nominal amplitude 

j He =0.3lb/sec=0. 136kg/sec (Table 2) were used for the simulation. 

One can see from Eq. 2-5 the time-trace pressure deviates from the telemetric data during 
the pre-pressurization stage. This divergence is a consequence of (i) the deviation of the actual 
pulse shape (Fig. 2-6) from the rectangular shape and (ii) the deviation of the experimental 
threshold pressure p th from its nominal value 20.54 psig used in the simulation. We 
incorporated the actual (experimental) rather than nominal value of the threshold and the actual 
pulse shape in the MNM simulation which shows an excellent agreement with the telemetric 
data (Fig. 2-7). In more detail, the pulses have amplitude j He =0.3lb/sec=0.136kg/sec and their 
form is given by: 

f l-exp(-(/-/,)/r ), t<t,+T 

m=jJ w (2D 

[exp(-(/ -T p )/ T lall ), />/, + T p 

where at t=t, the ullage pressure drops to the threshold value of p m and i-pulse is generated. 
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Table 2. Tank Ullage Pressure Dispersions 


PARAMETER 

NOMINAL VALUE 

3a VARIATION 1 

Trajectory 

Normal & Early RTLS 


SSME 

#1 

#2 

#3 


• Pressure r 

3560 

3584 

3579 

+42.4/ -55.9 psi 

• Temperature t 

830 

820 

852 

+31 .9/ -39.9 *R 

♦ Orifice Area 

0.02871 

0.02761 

0.02944 

±2% 

♦ Nominal Power Level 

104.5% 

104.5% 

104.5% 


• L0 2 Flowrate 

933.28 Ibjsec 

+6.33/ -6.78 Ibjsec 

Ullage Volume at End of 
Replenish 

137ft 3 

±22.3 ft 3 

Pre-Pressurization 





• Transducer Accuracy 

Tavis Flight Experience 

±0.4 psi 

• Drain Back / Seal Leak 


21.4 Ibjsec 


± 1 .5 lb Jsec 

* Max. Hold Time 
He Pre-press Supply 


0 sec 


5.0 min. 

- Flowrate 


0.3 Ibjsec 


± 0.05 Ibjsec 

- Supply Temperature 


572.4’R 


+ 8.7/ -19.5 °R 

- Valve Open Time 

! 

0 5 sec 


±0.1 sec 

Ullage Pressure at T 9 

20.54 psid 

+ 0.176/ -0.348 psi 

Ambient Temperature 

79.75°F 

48.8/ 81.6 T 

L0 2 Splash Effects 

• Splash Coefficient 

4.375* 

±0.35 


Gas Mixture. 


Wall 


Temperature 


I 


316.91 

302.73 

288.55 

274.36 
260.18 

246.00 
231.82 
217.64 
203.45 
189.27 
175.09 

160.91 

146.73 

132.55 

118.36 
104.18 

90.00 



Gas -LOX Interface 


[K] 


Temperature 



111.69 

110.34 

108.98 

107.63 

106.27 

104.91 

103.56 

102.20 

100.85 

99.49 

98.13 

96.78 

95.42 

94.07 

92.71 

91 36 

90.00 


IK] 


Fig. 2-3: Stratification of the ullage and wall temperature of GOx/He during the pre- 
pressurization. Results of 3D simulation. 
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Y[m] 


Fig. 2-4: Temperature distribution along line 1 on a tank section at tpl pKSS = -138sec 
corresponding to the pressure p=2.42atm (see Fig. 2-5). Results of the 3D simulation. 



Fig. 2-5: Ullage pressure during the pre-pressurization and the first part of the flight: 
STS115DB telemetry (red) and 3D simulation (blue). 



Fig. 2-6: The form of the He mass flow pulses. 

We chose a single fitting parameter of MNM, numerical coefficients /3 U in the velocity 
correlation for the turbulent natural convection (see Eq. (22)), in such a way as to fit the 
temperature distribution found in the 3D simulation (Figs. 2-3 and 2-4). The sawtooth 
modulations of the ullage pressure induced by the He pulses were fitted to the telemetric data. 
For the fitting we chose parameters of He impulses in Eg. (21) (the pulse duration t and its 

pulse 

front and tail durations t and t ). 

front tail 


16 


h'; r =0.15 P h K UG oc (3 h x\ Ra x =- _ T ; , )_£ , 


f^L.G^L.g 


h\x = P u u N ( x ), u A x i) = 0-178 

x. 


Ra„ Pr 


NV2 


1 + 0.494 Pr 273 


°c P u x oru N (x i ) = 


SPrfLn 
. Cl,gPi.„ 


,1/3 


( 22 ) 


s J 


j N = 7nR^5 xUf4x p U G cc p u x im \ 

8 x = 0.565* [(l + 0.494 Pr 2/3 ) / Ra x Pr 23/ ’ 5 ] oc x V10 

The fitting gave the following values for the parameters: 

• |3 U = 0.05 

• t = 0.09sec, t = 0.5sec, and r . ~0.7sec. 

front tail P“ lse 

• The threshold pressure for a He pulse generation p th is chosen to fit experimental data 
shown in Fig. 2-7. 

Result of fitting the ullage pressure is shown in Fig. 2-7. We see that results of the MNM 
simulations fit very well the STS-115 data. We reemphasize that the only fitting parameter for 
the MNM, was the numerical coefficient (3 U in the correlation relation for the natural convection 
velocity. This coefficient depends on properties of the fluids and is essentially independent of 
the tank design. Therefore we used MNM with this fitting parameter for analysis of the LOx SLS 
tank. 
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° -150 -100 -50 0 

Time, sec 


Fig. 2-7: Ullage pressure during the pre-pressurization: STS115DB telemetry (red) and 
MNM simulation (blue). Results of the fitting with 6 U =0.05, t = 0.09sec, t = 0.5sec, and 

front tail 

t = 0.7sec. 

pulse 

Result of fitting of the MNM simulations of the ullage pressure teeth induced by the He 
pulses for other Space Shuttle flight is shown in Fig. 2-8. In order to fit the STS135 telemetry 
data He pulse parameters were changed only ~ 10% in comparison with the STS 115 values 
(Fig. 2-7). 
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Time, sec 

Fig. 2-8: Ullage pressure during the pre-press (-280sec<t<-160sec): STS135DB 
telemetry (green) and MNM simulation (blue). Results of the fitting with (3 U =0.05, T f o = 

0.09sec, t = 0.45sec, and t = 0.75sec. 

tail pulse 

The fitting parameters of the MNM are further used in modeling the pre-pressurization and 
fight regimes to assess the risks for the proposed design of the SLS tanks (Section 5). 

Numerical MNM and 3D simulations of the time-traces of the temperature of the ullage and 
wall at some regions of the tank are shown in Fig. 2-9 for (3 U = 0.05 and rectangular He pulses of 
nominal amplitude j He -0.3lb/sec=0.136kg/sec and duration 0.5sec. One can see that the 
corresponding curves for MNM and 3D simulations are closed enough (see also below Fig. 3-6). 
These results substantiate choose the value of p u = 0.05. 




Fig. 2-9: Time-traces of temperatures of ullage (upper line) and tank wall (lower line) in 
various layers in the tank: green curve corresponds to the layer near tank top. Results of 
the MNM simulation (left) and 3D simulation (right). 


3 Operation during the fight 


The GOx is used to cool the nozzles, where it is heated to the temperature T G ox= 
800R=444K. The hot gas is further injected into the LOx tank to keep the pressure at the 
required value which would otherwise decrease due to the LOx level lowering during the fight. 
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The rates of the LOx propellant and hot GOx flows controlling the ullage pressure are shown in 



Fig. 3-1: Flow rates of hot GOx into the tank and LOx from the tank during the flight. 


To fit the STS-115 data [2] the form of dips seen on the curves in Fig. 3-1 was slightly 
modified for the MNM simulations (~ 10%). The rates correspond to the nominal values 
presented in Fig. 3-1. Result of fitting the ullage pressure and temperature are shown in Figs. 3- 
2 and 3-3, respectively. The ullage pressure falls within (20.5-5-14.5) psig range and the 
temperature -- within (-100+1 00)F=(200+310)K range. We see that the results of MNM pressure 
simulations fit fairly well the STS-115 data (Fig. 3-2). 

The deviation of the simulated temperature from the data is substantial during the flight but 
small at the pre-pressurization stage (t<0, Fig. 3-3). In the latter case, as can be seen in Fig. 3-4, 
the width of the transition region for the temperature is very small (< 2cm) and we conjecture 
that the temperature sensor is located outside this region. The observed deviation of the 
experimental and simulated temperature during the flight (Fig. 3-3) can be attributed to location 
of the detector in the relatively cold transition region which is larger during the flight, ~ 5cm (Fig. 
3-4). Another possibility for the deviation is suggested by the 3D simulations (Fig. 3-5) which 
shows pocket of relatively cold gas in trapped in the upper corners of the ullage, where 
presumably the detector was located. These pockets are absent during the pre-pressurization 
(Fig. 2-3), which is a consequence of the relatively low density of He compared to GOx. 

Time traces of the ullage and wall temperatures at various locations in the tank are shown in 
Fig. 3-6. They correspond well to the pronounced temperature stratification detected in the 3D 
simulation (Figs. 3-5 and 3-6). The 3D simulation also shows that the initial complex 
temperature distribution (Fig. 3-5, t=2sec) is transformed into almost flat layers after lOsec. 
This result justifies the approach adopted in the MNM. 



-150 


-100 


-50 0 

Time, sec 


100 


100 


-50 0 

Time, sec 


Fig. 3-2: Time-traces of ullage pressure during the pre-press and the initial stage of 
the flight in relative (left) and absolute units (right). STS115DB telemetry (red) and MNM 

simulation (blue). 
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Fig. 3-3: Time-trace of ullage temperature during the pre-pressurization and the initial 
stage of the flight near the tank top where the thermo-sensor is located: STS1 1 5DB 
telemetry (red) and MNM simulation (blue). The green curve shows time-trace of the wall 

temperature. 



Fig. 3-4: Ullage temperature radial distribution at the level of the thermo-sensor 
location. Results of 3D simulation. 
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Fig. 3-5: Temperature and velocity distributions at t = 2sec and t = lOsec after the start 
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Fig. 3-6 shows that the increase of the wall temperature lags behind the growth of the ullage 
temperature. The lag is T gw «100sec which confirms theoretical predictions based on the 

estimates in Sec. 1.5. Time-traces of the partial pressures of GOx and He are shown in Fig. 3- 
7. 
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Fig. 3-6: Time-trace of temperatures of gas ullage (left) and tank wall (right) in various 
layers in the tank: red and blue curves correspond to the layers near tank top and gas- 
liquid interface, respectively. Results of the MNM simulation. 
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Time, sec 


Fig. 3-7. Time-traces of partial pressure of gaseous Ox and He in ullage layer near the 
liquid-gas interface. Results of the MNM simulation. 


Dynamics of the condensation rate and the condensed mass are shown in Fig. 3-8. 
Evaporation is seen to override condensation during the first stage of the flight (0<t<10sec, 
Fig. 3-8). This result is explained by the stratification of partial pressure of Ox and He gases and 
the dynamics of the ullage pressure. Indeed, the ullage pressure drops sharply during of the first 
stage of the flight (Fig. 3-2) following a rapid decrease of the liquid’s level. The partial GOx 
pressure drops and approaches latm (Fig. 3-7, left). As a consequence, the interface 
temperature T s decreases to the value T L (see Eq. (3)) and the condition Eq.(9) breaks down, 
implying that evaporation becomes dominant. Later on, the partial GOx pressure increases 
(Fig. 3-7, left) due to the GOx pumping. As a result, evaporation becomes overridden by 
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condensation (Fig. 3-8). Total mass of the condensed M cond ^ 1kg (Fig. 3-8, right), i.e., it is 
negligible compared to the total mass injected hot GOx: M GOx =1.3kg/sec x 100sec= 130kg. This 
means that in contrast to the pre-pressurization stage the contribution of condensation- 
evaporation process to the overall mass and heat balance during the flight (t>0) is very small. 


O 



Fig. 3-8. GOx condensation rate (left) and total condensed mass (right). Results of the 

MNM simulation. 

4 The limiting case of fully mixed gases 

The pronounced stratification of both ullage temperature and partial pressure of GOx and He 
was found in the 3D simulation (Fig. 3-5). The stratification is also observed in the MNM 
simulation provided one uses the value /3 U = 0.2 in the correlation for the velocity given by Eq. 
(22). This value of p u is chosen in the MNM simulation to best fit the average velocities given by 
the 3D simulation. It should be noted that the value of j8 u for small variations of the gas density 
where the correlations Eq. (22) strictly hold is /3 U ~ 1. For large temperature gradients near the 
wall the variations of the gas density may be substantial. Therefore, it is instructive to consider 
the ullage gas dynamics for (3 U = 1, where due to higher velocity of natural convection near the 
wall, the temperature stratification almost totally vanishes (Fig. 4-1) and stratification of partial 
GOx and He pressure is very weak (Fig. 4-2). 



Fig.4-1 : Ullage temperature time-trace at various locations in the tank: red and blue 
curves correspond to the layers near tank top and interface, respectively. Results of the 

MNM simulation for /3 U = 1 . 
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Fig. 4-2: Time-traces of partial pressure of GOx and He in the ullage; red and blue 
curves correspond to layers near the top and the bottom of the tank, respectively. 
Results of the MNM simulation for /3 0 = 1. 

However, even in this limiting case of total mixing of the two gases, the time traces of the 
ullage pressure and temperature do not change considerably (see Fig.4-3 ^-Fig. 4-5). 
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Fig. 4-3: Time-traces of ullage pressure. STS115DB telemetry (red) and MNM 
simulation (blue). Results of the MNM simulation for )3 U = 1. 



Fig. 4-4: Time-traces of ullage pressure; blue for stratification (/3 U = 0.2) and green for 
total mixing case (/3 U = 1). Results of the MNM simulation. 
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Fig. 4-5: Time-trace of ullage temperature near the tank top: blue for stratification 
(P u = 0.2 ) and green for total mixing case (/?„= 1 ). Results of the MNM simulation. 


In the case of total mixing the average condensation rate and, accordingly, the total 
condensed mass are smaller than those for stratification (Fig. 3-7, j 3 U = 0.2) during the late pre- 
pressurization stage (Fig. 4-6, t<0). On the other hand, they are essentially identical during the 
early pre-pressurization stage and the flight (Fig. 4-6, t>0). The decrease of the condensation in 
the absence of stratification during the pre-pressurization stage is explained by the increase of 
the gas temperature and decrease of the partial pressure of GOx near the interface compared 
to dynamics with stratification. The partial pressure decrease leads to decrease of the surface 
temperature T s in the right-hand-side of inequality(9), while the rise of the gas temperature 

T g makes the left-hand side of the inequality larger. Both lead to inhibition of condensation and 
acceleration of evaporation, which explains the observed dynamics (Fig. 4-6). 



Fig.4-6. GOx condensation rate (left) and condensed mass (right). Results of the 

MNM simulation for /3 U = 1. 
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5 Scaling of the risks 

5.1. Scaling of pressure near the feedline top and bottom 


The pressure in LOx tank feedline is high due to hydraulic pressure of a large mass of LOx. 
The total pressure near feedline top p top and bottom p bo t are 

Phonon, = SPl ( ,a„k “ L «Uage + L p,pe ) + Put, age =5.9 atm (t < ~\ 5(k, P ullage = P am ) 

/ \ (23) 

Ptop = SPl (^lank - ^ullage ) + P ullage = 2 - 56 atm (5.45a/m at max acceleration), 

The pressures during pre-pressurization for 3 and 5 engines are approximately equal to 
(2.6(3.6)atm at t < -152s, p uUage * latm 

1 4(5 )atm at - 1 38s < t < 0, p ullage « 2.45atm ’ 

5.9(92)atm [t < -152s, p ullage « latm) 

K, < 

7 3(10.6) atm (-138s <t <0, p ullage «2.45atm) 


top 


Pbo 


(24) 


The Shuttle achieves maximal acceleration grocket := 2.5go at the point t=80sec, when the liquid 
level is 1.69m above the tank bottom (Fig. 5-1). The maximal top and bottom pressures 
corresponding to total acceleration g = 3.5g 0 are approximately equal 
p lop =53(S.5)atm at t>50sec, p ullage »1.75atm, 

Pbouom =17.5(28 )atm (t>50sec, Pullage * 1.75atm) (25) 

For estimation of the hydrostatic pressure in the feedline for the SLS tank with 5 engines we 
assume that the tank height and pipe length scale as 5/3Htank and 5/3L pi p e (H^k =1 5.05m and 
L P i P e =29.4m are the parameters of Shuttle LOx tank). 



Fig. 5-1 Time traces of Shuttle acceleration (left) and the liquid level in the tank (right). 

Thus, the maximal pressure at the bottom of the feedline for SLS with 5 engines will reach 
28atm. Such high pressure in the pipe and the engine's inlet has to be taken into account for the 
development of SLS. 


5.2. Scaling of shape and size of the LOx tank for SLS with 5 engines 


Two possible shapes of the SLS tanks are considered in this work: (i) the Shuttle LOx-tank 
shape and (ii) the Shuttle LH2-tank shape with volume upscaled to 5/3 of the Shuttle LOx tank 
(Fig. 5-2). The diffuser size is changed in the case (ii) (Fig. 5-2) to prevent contact with the liquid 
surface. The initial ullage volume V 0 and the volume in end of pre-pressurization i.e., at the 
points t=-152sec and t=-1 38 sec (see Fig. 2-1), respectively, are increased by the factor 5/3: 
V 0 =8.43m 3 and 13.2m 3 . The fluxes of hot He and GOx are also increased by the factor of 
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5/3 compared to the Shuttle. Fig. 5-2 also shows the gas velocity field near the diffusers for both 


shapes of LOx S 


S tanks. 



Fig. 5-2 The shapes of SLS tanks and diffuser considered in this work: (i) the Shuttle- 
LOx-tank shape (left); (ii) the Shuttle-LH2-tank shape (right). 


Predicted time-traces of the ullage pressure and temperature for the two shapes of SLS 
tanks are shown in Fig. 5-3 and compared to the Shuttle nominal time-traces. 



Fig. 5-3 Predicted time traces of ullage pressure and temperature for different shapes 
of the SLS LOx tanks. Results of MNM and 3D simulations. 
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Fig. 5-4 Predicted time traces of ullage pressure for different shapes of the SLS LOx 
tanks. Results of the MNM simulation. 



Time, sec 


Fig. 5-5 Predicted time traces of ullage pressure for LH2 shape of the SLS LOx tanks 
with different initial ullage volumes. Results of the MNM simulations. 

One can see that the time-traces of ullage pressure for the two LOx tank designs are close 
to one another. It is important to note that the predicted pressure time-traces lie within the band 
of acceptable ullage pressure (Fig. 5-4). However, the pressure time-trace for the LOx tank of 
shape (ii) lies very close to lower bound of the band, creating hazardous situation. To analyze 
the possibilities of controlling the risks we carried out the simulation for several other initial 
ullage volumes and showed that the situation can be substantially improved if the initial ullage 
volume in the LOx tank is decreased by 10% (Fig. 5-5). 

6 Effect of sloshing in the LOX tank 

Intensive sloshing of cryogenic liquid (LOx) can occur after the start due to vibration of the 
rocket. Sloshing of cryogenic liquids can change the ullage pressure due to variations of the 
heat and mass fluxes at the liquid-vapor interface and near the tank wall. Two main types of 
sloshing are shown in Fig. 6-1. Generally, the motion of the LOX tank during the lift-off can be 
decomposed into three major components: twang, breathing, and oscillation. The twang occurs 
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before the lift-off, while the breathing and the oscillation occur during the lift-off. The twang 
occurs approximately 5 seconds before the lift-off when the main engines are ignited and the 
Shuttle is released from the launch pad. When the main engines are ignited, a torque is created 
around the attachment point of the rocket boosters which makes the entire vehicle tilt [9]. The tip 
of the ET tilts as far as 30 inches from the rest position and then returns back. When the tank 
returns to the near vertical position the solid rocket boosters are ignited and the launch takes 
place, 


Fig. 6-1 Two main types of sloshing: (i) Longitudinal sloshing waves excited by lateral 
forces (left); (ii) Axial sloshing (capillary waves) excited parametrically by axial vibrations 

(right) 


This motion initiates oscillations of the liquid surface with period T osc ~2.5 sec and amplitude 
a osc <0.5m (Fig. 6-2). 



Fig. 6-2 Tilting of the Shuttle from the vertical direction after the lift-off (left) and 
typical sloshing oscillations (right). 

6.1. Longitudinal sloshing waves 


The rocket oscillation leads to variation of the liquid level near the tank wall [12,13]. Our 3D 
simulations showed that the wall temperature near the liquid-vapor interface is lower than the 
boiling temperature (T b0 ii>1 00K at p>2.25atm) for distance from the interface less than 0.75m 
during both pre-pressurization and the flight (see Fig. 6-3). It means that the variation of the 
liquid level near the tank wall will not stimulate evaporation of the liquid provided the oscillation 
amplitude a 0S c is smaller than 0.75m. This condition is apparently satisfied even for the Shuttle 
(see Fig. 6-2). 

The SLS design differs from the Space Shuttle construction in that no torque on the rocket is 
created by the engines ignition. As a consequence the hazard associated with the twang of the 
ET is essentially eliminated in the SLS design. 


28 



180 


^ 160 
0) 

3 

03 


E 

0) 

TO 

5 120 


100 


Distance from interface along the wall, m 


Fig. 6-3 Distribution of the wall temperature near the liquid-vapor interface for different 
times (x=0 is the interface position for a given time). Results of 3D simulation. 


After the start (t > 0) the stationary longitudinal sloshing oscillating waves enhance 
condensation due to (ii) increasing the area of the liquid surface and (mainly) due to (ii) mixing 
the heated surface liquid layer with adjacent cold LOx with T=T L . As a result of mixing the 

heated surface layer became narrower but the interface temperature did not change, it is given 
by the partial GOx pressure: T = T s (p GOx ) =100K. The condensation flow rate at the liquid-gas 

interface is limited by heat balance between the hot gas T g near the surface and the liquid 
surface with T=T : 


J com! ) 
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1 
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1 

* L L (t) ^ L g (t ) 

1 

+ 


(26) 


The maximum value of the condensed mass M cond induced by the intensive sloshing can be 

estimated assuming complete mixing of the surface layers of cold and warm liquid over each 
period of oscillation. Due to such intensive mixing the width of the heated liquid layer over a 

fraction S',* of the surface area will be very small at end of each oscillation, L l — >0 (Fig. 6-4). 
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Fig. 6-4. Sketch of the temperature distribution before and during the sloshing 


In this case the first term in the right part of Eq. (26) is much larger than the second. 
Therefore, the mass condensed for one period can be estimated as 
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(27) 


= f «S;(T, -T,) 


4K,C lPl T mc 
Ql + C H 2 {j g ~T^s) 



The total condensed mass induced by the longitudinal waves over the time of t is given by 


< 2(0 = KoMJ 


As,+,sy 




T , 

V 1 OSC J 


The total condensation flux is equal to 

vJ sbsh _ max = q 26 i£ 


maXj cond 


sec 


(28) 


(29) 


forS* * S L / 3, T(p GH2 ) = 100 A:, 100 AT <T g < 200 K near the interface 

It should be compared to the much larger the flux of the injected hot GOx: 

J go* = 2 ■' 6 ik x [ N nozzle 7 3 L N nozzle ~ number of nozzles in SLS 

Thus, the longitudinal sloshing waves cannot appreciably change the ullage pressure. 

Effect of sloshing is expected to be even smaller for the SLS since the flow rate of injected 
hot GOx will increase with larger nozzle number. However, SLS vehicle has higher thrust and 
one can expect stronger vibrations due to adding one or two RS-25 engines. 


6.2. Axially excited sloshing: generation of droplets 


The thrust oscillations initiate axial excitation and lead to Faraday instability of the liquid 
surface. Intensive oscillations of these modes can result in the liquid level discontinuity and 
generation of the liquid droplets under the liquid-gas interface. This process will change heat 
balance on the interface. Evaporation of the cold liquid droplets into the hot ullage gas can 
increase the ullage pressure during the beginning stage of the flight. At the same time enlarging 
the cold surface can accelerate condensation. Since the rate of evaporation of liquid Oxygen 
depends on partial pressure of gaseous Ox, in the case where He is used as the pressurization 
gas we expect enhanced evaporation and consecutive increase in a ullage gas pressure. For 
gaseous oxygen pressurization these droplets will enhance the condensation rate. 

Dispersion relation for small wavelength capillary wave can be presented as [12,13]: 

f „ \ 


co = 


gk- 


P J 


tanh(A/7) 


where a is surface tension, k=2TT/A, A as the wavelength, k is the wave number and p is the 
liquid density. At small k«1 

,2 _g^ BL tanh^ft/a) (30) 


co. 


a 


The nonlinear dynamics of the droplet formation and ejection is very completed and strongly 
depends on the rocket vibrations spectrum. Here we carry out a simplified analysis this 
dangerous effect. The typical diameter of ejected droplets is about [14,15] 

djrop - 1' 7 ! / P L ) m co' m , 0) = 2 71 f (31) 

where p L =1141kg/m 3 and o L =0.0135N/m are liquid density and surface tension, co is 
vibration frequency (co =2xf vib r). Characteristic values of d drop are shown in Fig. 6-5 


(left): d u (2.5Hz) - 0.8 mm and c/ (200 Hz) — 0.2 mm . 
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Fig.6-5. Typical diameter of ejected droplets depending on thrust oscillation 
frequency (left) and generation of droplets under action of supercritical vibration 

acceleration E d >0 [13]. 


Each projectile area of diameter d pr0 j ~5d drop ejects the droplet with a rate given by [14,15] 


dn/dts h = 0.04 s d \o [s~'],e d = (a vibr -a d )/a d . 


1/3 ,.4/3 


a d = 0 . 26 ( 0 -/ pTg) 


(32) 


Here a, co are vibration acceleration and frequency, a d is the critical driving acceleration for 
the onset of droplet ejection. The droplet ejection occurs if vibration acceleration a vjb r >ad, i.e. 
when the frequency and amplitude vibration A v ibr satisfy the condition: 


«v,v = Aibr™ > a d or >4=0. \mm(65Hz / f jbr ) 


213 


(33) 


4 = 0.2mm for f = 25 Hz, A c = 50 jum for f = 20 0Hz 

Sufficiently strong vibrations can occur only after the start. At this point the thermo-diffusion 
length in the liquid has increased to the value of l_ D (100sec)=7mm during the prepress (t = 

lOOsec). This value is much greater than the typical drop 
diameters: L D — 5mm » d Jnip for f > 50Hz) . It follows that inner temperature of generated 
droplets T drop is close to the surface temperature of the liquid T s (p GOx ). Therefore, the liquid 

mass evaporated due to the vibration is very small. Indeed, the mass evaporated with surface 
of single droplet s during the droplet life-time t drtfe can be estimated as 
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f/j 

evap^d drlife) ^ drop 


pros + 1 drlife 
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J cond (tpreJS drop 


’ ^ drlife * 


2 h 


drop 


1 drlife 5 1 drlife 
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S l ' ' V S 

where h drop is the maximum height an ejected droplet rises above the surface (see Fig. 6-5). 
The total evaporated mass is given by 
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Here dN'ff / dt is total number of droplets ejected over unit time from the liquid surface S L 


drop 


Taking into account Eq. (34) and condition d 


pro] 


>5d dkvi s ltv ls ml =AI2S we conclude from 


proj 


Eq.(35) that the total evaporation mass is a negligible quantity 
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5 ) <0.03g/sec 


evap \ prepress > 

In the estimation we used the results of calculation presented in Sec. 2 (see Fig. 2-8) 

^(^„)^0.2g/sec. 

6.3. Combined action of longitudinal and axial sloshing vibrations 

An appreciable effect can occur only as a result of combined action of intensive longitudinal 
and axial sloshing vibrations (Fig. 6-4, left)). The maximum value of the condensed mass 
induced by sloshing waves can be estimated assuming intensive mixing of surface layers of 
cold and warm liquid during each period. Due to such mixing the thickness of some area of the 
heated surface layer S] will very small at end of each oscillation, L l -*0 (Fig. 6-4). Droplets 
generated by the axial sloshing will have T near T l , except for a very thin surface layer l_ L — >0 

near the surface having T=T s (p G ox)- In this case the maximum value of the condensed mass 
m on the surface s, of a single droplet during the droplet life We can be estimated as 

cond drop 
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cond \ drlife 


drop 
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The total mass flow rate of condensed vapor on all ejected droplets is 

j total _ ( \ dN drop _ f 4 S I m(t drlj f e )n d r tip 
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Using S, ~S L /3 and t drUfe 


rate of injected hot GOx mass flow 
^ GOx ~~ 2.6(AU,/3)i, -number of nozzles in SLS . 


drop i g * 0.02sec we obtain from Eq.(38) that the total mass 

flow rate of condensed vapor on all ejected droplets is < 0.45£g/sec . We notice that the 

higher: 


is 


much 


As a result of the condensation the ejected droplets will be heated and falling back into the 
liquid will increase the temperature of the liquid surface layer. The width of the heated surface 
liquid layer can be found from the heat balance equation 

45 [il * c oo.{t s - r,)]-» = c, Pl (t s -t l )s l - C2, 


' _ ''t \t i 1 : ■ f _ fTT 

1 heated V / ^ ^ 25/7" j' drop* 1 dr life W ^ > 

mt)- 5L iM = 5cm for 'w * — > W * °- 02sec 


( 39 ) 
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Thus, combined action of longitudinal and axial sloshing vibrations can lead to formation of 
relatively thick hazard surface layer of heated liquid with temperature T-23.8K, thickness 

L*5cm, volume, V hea =S t3Jxk I%? aed l 2*1 Am 1 2 3 4 5 and mass M ka/ = V hut p L « 1600kg . Total mass of 
LOx in the tank is 627500kg. 

SLS vehicle has higher thrust and expected vibrations due to adding one or two RS-25 
engines. At the same time all the effects under consideration are taking place in a small layer 
adjacent to the interface and the mass flow rate of injected hot GOx increases with the engine 
number. Therefore the influence of these effects on both the mass and the heat balance in the 
tank will not be critical. However, the effect of heating of relatively thin layer to T«100.8K can 
arise in LOx SLS tank. 

SLS vehicle has higher thrust and expected vibrations due to adding one or two RS-25 
engines. At the same time all the effects under consideration are taking place in a small layer 
adjacent to the interface and the mass flow rate of injected hot GOx increases with the engine 
number. Therefore the influence of these effects on both the mass and the heat balance in the 
tank will not be critical. However, the effect of heating of relatively thin layer to T=100.8K can 
arise in LOx SLS tank. 


Conclusion 

We performed theoretical investigation of main physics processes in LOX tank during pre- 
starting prepress and rocket flight that can lead to possible hazards in the case of rescaling of 
main parameters of Shuttle with regard to proposed LOX tank designs for SLS with 5 engines 
(the situation with 4 engines is less critical). Our investigations based both on theoretical 
analysis and on MNM and 3D ANSYS package simulations allow us to make the following 
conclusions. 

1. Hazard analysis for LOx ET shows that the pressure control bandwidth remains almost 
the same for the SLS. The time-traces of SLS LOx tank ullage pressure and temperature were 
calculated after an appropriate rescaling of the relevant parameters. The rescaling takes into 
account the increase by factor 5/3 of the linear dimension of the tank compared to the tank of 
the Shuttle. The rescaled parameters include flows of hot He and GOx, and also initial volumes 
of ullage and liquid in the tank. The predicted time-traces of the ullage pressure are found to fall 
within the acceptance limits for both SLS LOx tank shapes. 

2. The MNM can reproduce accurately the pressure during pre-press and blow down. 

3. For the GHe and GOx flow rates and the initial ullage volume rescaled by 5/3 for the 
SLS LOX ET, no new hazards are found for the tank performance; the temperature stratification 
is more pronounced in the rescaled SLS tank. Such stratification occurs due to slow mixing of 
GOx/He gases: the heavy cold GOx localizes near the liquid surface and the injected hot light 
He accumulates near the tank top. The stratification is predicted theoretically and confirmed 
and scaled by the MNM and 3D simulations. 

4. Parameters of the MNM, developed in our group, were fitted to make the MNM 
reproduce the ullage temperature stratification observed in the 3D simulations. Remarkably, it 
is found that an accurate correspondence between the model and the 3D simulations is 
achieved by fitting just a single parameter: the scale of the velocity of natural convection in the 
turbulent boundary layer near the tank wall. 

5. The MNM reproduces with high accuracy the telemetric data of saw-tooth modulations of 
the ullage pressure during pre-pressurization stage, which result from injections of He pulses. 
This accuracy only required fitting the shape of the He pulses. Likewise, the model is found to 
reproduce fairly well the telemetric Shuttle data of the ullage pressure variation during the flight. 
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6. The ullage pressure in LOx tank is shown to be essentially insensitive to temperature 
stratification: a model calculation implies that even in the case of complete mixing of the cold 
and hot gases, where the temperature stratification is absent, the ullage pressure time-trace 
does not change significantly. 

7. In contrast to the SLS LOx-shaped tank, the time-trace of the ullage pressure in SLS 
tank with LH2 shape is found to approach the lower acceptance limit. It is shown that this 
hazard can be eliminated by a decrease of the initial ullage volume in the tank. 

8. Theoretical analysis of bubbles dynamics in LOx has been analyzed. It was shown that 
generated bubble volume increases as the tank height squared. Collapse of the bubbles during 
the early stage of the pre-pressurization has been demonstrated. 

9. Condensation-evaporation flows are found to weakly affect the ullage pressure. First, 
the condensation is controlled by the high pressure of GOx in the ullage. In contrast to one’s 
naive expectations, even in a hypothetical case where the gas temperature in the vicinity of the 
liquid-vapor interface reaches the temperature of the injected He (320K) or GOx (800R=444K), 
it is condensation rather than evaporation which is dominant during the early pre-pressurization 
stage and the flight. In reality, due to stratification the actual temperature near the liquid-gas 
interface is much lower, which leads to strong suppression of evaporation during all stages of 
the pre-pressurization and the flight. Second, the condensation is controlled by the slow heat 
transfer rate in the LOx. It is predicted that the total condensed mass during the pre- 
pressurization and the flight is negligible. Therefore, large change of ullage pressure relating to 
condensation-evaporation process at the liquid-vapor interface is unlikely. 

10. The estimated temperature of the interface is lower than 100K and the estimated 
thickness of the heated liquid layer is smaller than 1cm. This layer corresponds to the total 
mass of the heated liquid which is less than 450kg during the fight (t fl ig h t =500sec). This mass is 
negligible compared to the total LOx mass in the filled tank, which is 627,500kg. Therefore, it 
does not play the role in the ullage temperature stratification and pressure dynamics. It was 
also found that the total heat associated with the condensation of GOx is much less than the 
heat of hot GOx injected during the flight. 

11. Slight temperature stratification of LOx is formation of thin warm liquid layers with 90< T 
< 98K near the interface and a layer with temperature T< 90.8K due to the liquid rising flow 
inside a very thin layer of thickness about 0.5mm. This liquid stratification does not play the role 
in the ullage temperature stratification and pressure dynamics and does not lead to breach of 
condition NPSP= P ul ,+P/ Kad —p Jhc —p svap (T s )>30psi for the net positive suction pressure due 

to the high hydrostatic head pressure p head near the feedline bottom. 

12. Intensive longitudinal sloshing waves induced by the laterally excited rocket vibration 
change weakly the temperature liquid and ullage stratification and the ullage pressure during the 
flight. Maximum rate of condensation flow arising due to the strong mixing of the warm surface 

liquid layer with the cold deep layers is limited by value of = 0.26^-. The mass flow rate 

of injected hot GOx is much more J GOx (t) = 2.6-!^ x(n noale / 3) , n nozz/e is nozzle number in the 
SLS. 

13. The vertical vibrations with the amplitude A vibr and the frequency f vibr such as 
Aibr > A = 0.1ww(100//z//j, lV ) 2 ' 3 induce generation of the warm droplets under the liquid-gas 

interface. Effect of axially excited vibration sloshing on LOx tank is extremely small as the 
temperature of generated warm droplets is very close to the interface temperature and 
condensation is depressed. 

14. Combined action of intensive longitudinal and axial sloshing vibrations can induce 
relatively great condensation and formation of relatively thick hazard surface layer of heated 
liquid with temperature T=T s (p GOx ). The rate of the condensation can 

reach J%£ = 0.45A:g /sec and thickness, volume, and mass of surface heated hazard layer can 
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reach L L =5cm, V=1.4m 3 , and 1570kg, relatively. Total mass of LOX in the tank is 
627,500kg x^/ 3). 

15. A possible hazard is related to the increased length of the tank and the feed line joining 
LOx tank with the engines in the SLS. The height of the tank for SLS with 5 engines is about 
25m and length of the feed line is about 50m. The hydrostatic head pressure near the bottoms 
of the tank and the feedline are 2.6amt and 9.2atm, respectively. The rocket acceleration during 
the fight reaches 2.5g 0 causing the pressures to increase to 17atm and 29 atm, respectively. 
Such high hydrostatic pressures have to be taken into account in the process of manufacturing 
the feedline and the engine pumps with LOx propellant. 
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